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Hadronic transitions from the lattice 

UKQCD Collaboration: C. Michael 8 

a Theoretical Physics Division, Dept. of Mathematical Sciences, 
University of Liverpool, Liverpool L69 3BX, UK 

I discuss strategies to determine hadronic decay couplings from lattice studies. As a check of the methods, I 
explore the decay of a vector meson to two pseudoscalar mesons with Nf = 2 flavours of sea quark. Although we 
are working with quark masses that do not allow a physical decay, I show how the transition rate can be evaluated 
from the amplitude for p — > nn and from the annihilation component of nn — > nn. I explore the decay amplitude 
for two different pion momenta and find consistent results. The coupling strength found is in agreement with 
experiment. I also find evidence for a shift in the p mass caused by mixing with two pion states. I also present 
results for the decay of a hybrid meson, for the case of heavy valence quarks. 



1. INTRODUCTION 

Rather few hadronic states are actually sta- 
ble under strong interactions. Because of this, 
a full understanding of hadrons from QCD will 
involve developing techniques to treat unstable 
hadrons and to evaluate their hadronic transition 
strengths. 

An area of hadronic physics that is of consid- 
erable interest is that of gluonic excitations: the 
exploration of hadrons with a non-trivial gluonic 
content. The most clear case is that of glueballs, 
however, any glueball states will mix with quark- 
antiquark mesons in practice. This mixing can be 
thought of as mediated by a hadronic transition. 
Moreover, the eventual hadronic decays of such 
mixed mesons also need to be evaluated to give 
a clearer picture of the relationship between the 
physical spectrum and the ingredients (glueballs 
and quark-antiquark mesons). 

The other natural area to explore gluonic exci- 
tations is that of hybrid mesons. These are quark- 
antiquark mesons with an additional non-trivial 
gluonic component. The benchmark example is 
a spin-exotic meson (so-called because it cannot 
be made from combining a quark and antiquark). 
These spin-exotic hybrid mesons will also be un- 
stable, and study of their decay transitions will 
guide the interpretation of experimental candi- 
dates. 

An understanding of methods to treat purely 



hadronic transitions will also help in developing 
methods to explore weak decays to hadronic final 
states (e.g. K — > 7T7r). 

The lattice QCD approach is the only quanti- 
tative non-perturbative method but the formula- 
tion in euclidian time causes problems with the 
treatment of decays as I now discuss. 

1.1. Euclidian time 

Hadron masses are extracted from lattice QCD 
calculations using two point correlators, however, 
their behaviour as e~ Mt in euclidian time will 
not be appropriate for hadrons that can decay. 
The most naive modification of the lattice QCD 
formalism caused by the introduction of decay 
widths is the replacement M — > M — iT/2 where 
P is the decay width. This gives a behaviour as 
e -Mt e iri/2 ^ mg j g na (. cong i g +;ent: the correla- 
tion must be positive definite yet this expression 
oscillates. Here I discuss this apparent paradox 
more fully pQ. 

The euclidian time propagator of an unstable 
particle which can decay into two stable particles 
of mass fj, can be expressed as 



1 r°° 

G(t) = - dE 

K J 2)i 



-Et 



P{E) 



(1) 



where we can express the spectral density for a 
resonance of mass M as 
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P( E ) = 7TT7 



T{E)/2 



2M (M - E) 2 + {T{E)/2f 



(2) 
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As an illustration, consider taking T(E) as a 
constant with T < (M-2/z) and (M-2p)t > 1. 
Then we have, in this approximation, 

e-2/itp 



2MG(t) 



-Mt 



cos(R/2) 



(3) 



2tt(M - 2^i)H 

This expression, when evaluated, shows that the 
threshold contribution (the second term) domi- 
nates at larger t and hence the oscillating be- 
haviour of the first term will be obliterated. This 
clarifies the apparent paradox mentioned above 
concerning the sign of G(t) : the contribution from 
the lower lying (in energy) component of the res- 
onance is enhanced in euclidian time. Indeed this 
can be seen directly from eq. ^ 

By contrast, the Minkowski behaviour of the 
two particle correlator under the same assump- 
tions is 



2MG M (t) 



= p -iMt -Tt/2 



ifj,t~p 



(4) 



2ir(M - 2p)H 

for which the second term is almost always negli- 
gible. 

In principle, a lattice determination of the spec- 
tral function p(E) from measurements of G(t) us- 
ing eq. would be sufficient to give a full de- 
scription of hadronic decays. This is not a simple 
task: the inverse Laplace transform need to ob- 
tain p(E) from G(t) is not numerically stable. It 
has to be stabilised by making model assumptions 
about p(E), for example as discussed in refs. jll2j . 
Another possibility is that the maximum entropy 
method may provide sufficient numerical stability 
and this has been tested in model cases 0. 

The same observation that the region near to 
threshold dominates is at the root of the conclu- 
sions [3] of Maiani and Testa that two body states 
have unappetizing properties in euclidian time. 

1.2. Particles in a box 

In practice this issue is less relevant since the 
lattice evaluations are performed in a finite spa- 
tial volume with periodic boundary conditions. 
This finite spatial size of the lattice implies that 
two-body states are actually discrete. By mea- 
suring their energy very precisely as the spatial 
volume is varied, it is possible to extract the 
scattering phase shifts and hence decay proper- 
ties. 



Consider the two body states of two pions of 
mass p with momentum ±k which will have to- 
tal energy E = 2(p 2 + k 2 ) 1 / 2 . Then if they 
are non-interacting, their momentum on a spa- 
tial hyper-torus of size L will be discrete with 
k = 2im/L where n has integer components (and 
for future use we define N such that there are N 
values of n having a given n = |n|, so N = 6 if 
n = 1). If the two pions interact with a finite 
range of interaction, compared to the spatial vol- 
ume, this implies that the finite size effect will 
result in an momentum shift of order L~ 3 7 pro- 
vided that L is large enough to contain one pion 
without undue distortion. The detailed analysis 
of Luscher |S] expresses the shift in this momen- 
tum value when there are interactions among the 
pions with phase-shift 8(k), provided E is below 
the inelastic threshold, as 



d{E 2 ) = 4d(fc 2 ) = 



-16nNtaxi5(k) 



(5) 



to leading order in L _1 , and indeed a much more 
precise expression with sub-leading terms is avail- 
able through Luscher 's work 

Now when there is a resonance at mn near to 
a two particle state, there will be an opportunity 
to deduce the resonance width. Then, one can 
describe the phase shift 8{k) by using the expres- 
sion for elastic irn scattering dominated by this 
resonance pole: 



tan<5(fc) 



T(k) 



2{m R - E) 



(6) 



so that determinations of S(k) at two or more 
values of k will enable the resonance parameters 
itir and T(kn) to be evaluated, where m 2 R = 
4(/i 2 + kft). The most promising way to deter- 
mine the two particle energy accurately is by us- 
ing several lattice operators (including some built 
like two pions and some built like the resonance 
under study) as sources and sinks to get the best 
determination of the energy levels from fitting 
the matrix of correlators. Note that in practice 
it will be difficult to obtain sufficiently accurate 
energy determinations of the lowest two-particle 
state and it will be even more difficult to deter- 
mine the energy level of the next heavier state 
with larger momenta. 
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1.3. A direct approach 

It should be possible to extract the hadronic 
transition amplitude directly from the lattice 
(rather than via energy level shifts), and this we 
now discuss. 

Using suitable lattice operators to create R at 
t = and annihilate a two-pion state with mo- 
menta k and —k at time t, the contribution to 
the correlator from a R state with mass m R and 
a two-body state with energy is given to lead- 
ing order in x by 



tl 



are-®* 



(7) 



where there is a summation over the intermediate 
i-value t\ and where h and b are the amplitudes 
to make each state from the lattice operators used 
and x is the required transition amplitude (R\tttt). 
Here we are assuming that the states R and tttt 
are normalised to unity. By obtaining h and b 
from the R — » i? and 7T7T — > 7T7T correlators, one 
can attempt to extract x. 

The complication, however, is that removal of 
excited state contributions is tricky, even if vari- 
ational methods are used to construct improved 
lattice operators to create the ground states. For 
example, if m R — E^ > then the transition 
time ti will be preferentially near (since the 
heavier state then propagates less far in time) 
and one can complete the sum over t% obtain- 
ing a ^-dependence of eq. 0as e~ £ " f . This same 
t-dependence would be obtained if the state with 
mass ttiji were to be replaced with an excited 
state with an even heavier mass. Thus one can- 
not separate the ground state and excited state 
contributions even in principle. See ref. |6I7) for 
a fuller discussion. 

The way forward is that if rriR = E^n, the 
ground state contributions to eq. have a t- 
dependence as te~ Ewwt whereas any excited state 
contributions behave as e~ E ** t as above. So now 
we do have a way to isolate the required ground 
state contribution: 



xt = 



Cfl-TTTr^) 



[CR-n{t) Ct, 



r(*)F 



-+0{x)+0(xH d ){8) 



Note that this separation is only by a power of 
t which is less than the case for two-point corre- 



lations where the excited state contributions are 
suppressed by an exponential e~( m «~" lfs '*. 

In practice the requirement of energy equality 
can be relaxed. Defining A = m R — E n7ri then 
the ground state contribution to the expression 
of eq. El evaluates to 2x sinh(tA/2)/(tA) = x(l + 
(£A) 2 /24 + ...). So this will be equivalent to the 
expression with A = provided 



(m R - E n7r )t < 5 



(9) 





Figure 1. Quark pair production (wiggly lines) 
for the three point function R — » n + n in euclid- 
ian time (running horizontally) . The left hand di- 
agram has the interpretation of a transition (our 
x) at an intermediate time while the right hand 
diagram can be thought of as some intrinsic mix- 
ing in the R state. 



So far we have described the behaviour of the 
Cn-^Trit) in the limit of small x. The first correc- 
tion term from multiple transitions (to eq. |S} will 
be of order x 3 t 3 /6 so we need xt <C 1. As well 
as these corrections which are of higher order in 
xt, one must also consider the intrinsic mixing 
of the initial R state with 7T7t (and vice versa). 
This intrinsic mixing (i.e. coming from the lat- 
tice operators used to create the states having 
an admixture - see fig. ^ f° r & n example) is ex- 
pected be of order x/£ where £ is the energy of 
the quark pair and so will contribute a term like 
xe~ Et /£. This is a contribution similar to that 
from excited states and so will be dominated at 
large t by the xte~ Et term we are looking for. So 
we need both xt to be small and t to be large. 
This implies that x must be small for this direct 
approach. And indeed, it is only when x is small 
that the more rigorous approach of determining 
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the energy shifts will be computationally difficult 
to implement. 

In summary, when t(mu — E vv ) <C 5; xt <C 1 
and (£" — E)t ^> 1 (with E' —E the energy gap to 
the first excited state), it is possible to estimate 
the hadronic transition amplitude directly, and 
this we will explore in detail for p decay. 

1.4. Comparison 

In order to show the relationship of this direct 
approach with Liischer's approach, consider the 
case when a two particle state is close in energy 
to the resonance. It should be possible to ar- 
range that energy levels are sufficiently close by 
choosing the lattice spatial length L. Then once 
x has been determined in this situation where 
the resonance does not decay significantly (be- 
cause the two particle level has the same energy as 
the resonance) , one can assume that the coupling 
strength determined can be used in a larger spa- 
tial volume where decay is energetically allowed. 

The most direct way to evaluate resonance 
decay is using first order perturbation theory 
(Fermi's Golden Rule) which implies a transition 
rate 

r = 2n{x 2 )p(E) (10) 

where the angle brackets indicate that an average 
over spatial directions will be needed. For a decay 
from the centre of mass with relative momentum 
k, the density of states p(E) = L 3 kE/(8n 2 ). 

Now to compare with Liischer's approach, we 
need to evaluate the energy shift of the two parti- 
cle state in a finite box (which is exactly N— fold 
degenerate) . This can be estimated from the mix- 
ing of the two energy levels which are close to each 
other, using second order perturbation theory (in 
x) which gives an energy shift, on the lattice, of 

Nx 2 

dE = -^-, (11) 

m R - E 

Thus when the resonance lies above the two par- 
ticle state, the mixing will move the resonance 
mass up and the two particle energy down, in the 
usual way as mixing causes repulsion of levels. 

Using the relationship between x and T from 
eq. we obtain 

4ttAT , . 

dE = - {m R -E)kEL> (12) 



which is exactly the same expression as would be 
obtained at leading order in L _1 from Liischer's 
formalism as presented above. This is as it should 
be, but does illustrate that an estimate of x di- 
rectly from the lattice may provide a useful way 
to evaluate these decay effects. 

Indeed the mass shifts are of second order in 
x so will be small in general, and hence difficult 
to determine with sufficient precision. The direct 
measurement of the transition giving x from the 
lattice is thus an attractive prospect. 

2. VECTOR MESON DECAY 

2.1. Transitions 

For on-shell transitions, it is possible to esti- 
mate hadronic transition strengths as described 
above and here we explore this approach further 
for the case of p meson decay to tttt, following 
ref. 

The situation we shall analyse is represented by 
the energy spectrum shown in fig.|3 here neglect- 
ing any interactions among the states. Note that 
no decay can actually take place with these lat- 
tice parameters. This is , however, optimum for 
our analysis since when the energy levels are ap- 
proximately degenerate, the hadronic transition 
amplitude can be evaluated more effectively. 

We evaluate correlations between lattice oper- 
ators creating both a p meson and a tttt state, 
using a stochastic method with sources on time- 
slices £ 9. to evaluate the quark diagrams shown 
in fig. |21 from 20 gauge configurations ^U] with 
Nf = 2 sea quarks of mass corresponding to 
about 2/3 of the strange-quark mass. These cor- 
relations, normalised by the two point functions 
as appropriate, are illustrated in fig. 0] Here the 
off-diagonal case (labelled p\ — ► shows the 

important feature that it grows approximately 
linearly with increasing t. As emphasized above, 
this linear growth will only occur for on-shell 
transitions and this is essentially the case here. 

To extract an estimate of the transition ampli- 
tude, we wish to evaluate x = {p\mr) where these 
states are normalised on the lattice (to unity) . If 
higher excited states are neglected, one can make 
a two state model (i.e. p and lightest irir state) 
with this transition amplitude and evaluate the 
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Figure 3. 



contributions as shown by the curves in fig. 
Because the location (at t\ between and t) of 
any transition is not known, one has potentially 
very serious contributions from excited states. As 
discussed above, these excited state contributions 
can be avoided if the transition is approximately 
on-shell, when the linear dependence on t of the 
off-diagonal transition (with value xt) is a unique 
signature of this on-shell transition that we wish 
to extract. 

Note that for the case of po — > ttitti, the on- 
shell condition is much less well satisfied but the 
relative momentum of the pions in the centre of 
mass is twice as large as for p\ — > t\qt\\ and hence 
x should be approximately twice as large (as we 
find), since for a P-wave decay there will be a 
momentum factor in the transition amplitude. 

A further check of the extraction of x comes 



from the box diagram, fig. 01, which will have a 
contribution behaving as x 2 t 2 /2 arising from a p 
intermediate state - see the data labelled ttittq — * 
7Ti7To in fig. 0| Again this contribution does seem 
to be present at the required level. 

Thus we have qualitative agreement from the 
three-point evaluations with two different relative 
momenta and from the four-point evaluation that 
xa ks 0.06. 




Figure 4. Ratios of 3 and 4-point correlators. 



A more quantitative estimate of x can be 
made 8 by reducing excited state contributions 
as much as possible. Thus one uses fits to remove 
excited state contributions from the normalising 
two-point functions and then evaluates the slope 
of the ratio of eq. |Sj versus t to remove the 0(x) 
term also. 

Note that the above methods can be used 
in the quenched approximation to estimate 
hadronic transition amplitudes. They depend on 
assuming that xa is fairly weak, as we indeed find. 

2.2. Energy shifts 

A more rigorous approach is to focus on en- 
ergy values. When two levels are close (our on- 
shell condition), then they will mix and the re- 
sultant energy shifts give relevant information [3] 
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as discussed above. Moreover we can estimate 
these energy shifts from our x-value which pro- 
vides more cross checks. These energy shifts can 
only be studied using dynamical fermions. 

From a full analysis (i.e. measuring all the 
quark diagrams illustrated in fig.[3]and fitting the 
matrix of correlators) we obtain the energy shift 
of the 7ri7ro state (i.e. the un-binding energy), 
as needed in Liischer's approach, as 0.02(2) (in 
lattice units) upward which is consistent but not 
sufficiently accurate to be of significant use. The 
energy shift of the p\ state can, however, be de- 
termined in this case because of a lattice artifact. 
The p with momentum 1 (in lattice units of 2it/L) 
can have its spin aligned parallel to the momen- 
tum axis (P) or perpendicular to it (A). Because 
the 7To7Ti state has relative momentum along a 
lattice axis and the transition from p to 7T7t has 
orbital angular momentum L=l (so a distribution 
like cob 9 in the centre of mass), only the parallel 
state (P) can mix with this two pion state. This 
situation is not typical - it occurs because the p 
with momentum 1 only has a transition to a tto'Ki 
state with a pion having momentum in the same 
direction. This mixing will not be present in the 
quenched approximation, so this provides a direct 
opportunity to see the effect of the two pion chan- 
nel on the p in unquenched studies. We do indeed 
find 8 a such a mass splitting between the P and A 
orientations of the p due to mixing for Nf — 2 but 
not for Nf = as shown in fig. |SJ Moreover the 
magnitude of this energy shift (0.026 ±7 in lattice 
units) is consistent with other determinations of 
the transition strength. 

2.3. Phenomenology 

The basic assumption is that the transition 
from p to 7T7T is given by an effective interaction 
with a finite spatial extent, this is usually sum- 
marised by an effective lagrangian where we nor- 
malise the coupling as g 2 = TME/k 3 in terms 
of the decay width. Then, provided the lattice 
spatial size is big enough that the hadrons are 
not distorted, our lattice situation (where no de- 
cay occurs) can be used to determine g and this 
can then be used to predict decay widths when 
the lattice volume is increased so that the mini- 
mum momentum becomes small. In our case we 
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will also need to reduce the quark mass to allow 
decay since we have m w /m p = 0.58. 

One complication is that the transition which 
is closest to being on-shell is p\ — > t^ittq which 
is not in the centre of mass. A generalisation of 
Liischer's approach has been made ^2 which al- 
lows for this. The more direct method we have de- 
scribed here can also be extended to this case jS] . 
There are possible problems in this treatment of 
decays of moving particles. For instance, on the 
lattice a boost to bring the p to rest will not bring 
the two pion state to have zero net momentum 
also. This is because the lattice situation allows 
energy non-conservation which, upon boosting, 
will imply some momentum non-conservation. 

From our lattice studies, we deduce that ax — 
0.06ti for pi — » ttittq. Translating [B] this lat- 
tice transition amplitude to the continuum nor- 
malisation, gives g — 

1.401 23- Usin S the ob- 
served pi energy shift gives another estimate, 
namely g = 1.561^3- These two values agree well, 
which is good support for our scheme of extract- 
ing hadronic transition strengths. 

Note that our lattice values would need to be 
extrapolated to light sea-quarks and to the con- 
tinuum limit to allow all sources of systematic er- 
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ror to be explored in the comparison with the ex- 
perimental situation. Nevertheless, these two val- 
ues agree well with the values extracted from ex- 
perimental data on decays of p, K* and <f> mesons, 
namely g w 1.5. This further underlies the via- 
bility of our methods. 

3. HYBRID MESON DECAY 

Hybrid mesons are those with non-trivial glu- 
onic excitations. The spectrum of such mesons 
has been determined from lattice studies, both 
for the case of heavy quarks and light quarks. 
Hybrid mesons have allowed couplings to meson 
meson systems. If these coupling are large, then 
the hybrid mesons decay widths will be big and 
such mesons will be difficult to identify experi- 
mentally. Here we present recent lattice results 
on these hadronic transitions. 

3.1. Hybrid mesons with heavy quarks 

One well defined way to treat heavy quarks 
(such as the 6-quark) on a lattice is to use the 
leading order of the HQET which corresponds to 
the static approximation. Here we take the heavy 
quarks to be fixed at locations R apart in the z- 
direction. 

Consider QQ states with static quarks in which 
the gluonic contribution may be excited. We clas- 
sify the gluonic fields according to the symmetries 
of the system. This discussion is very similar to 
the description of electron wave functions in di- 
atomic molecules. The symmetries are (i) rota- 
tion around the separation axis z with represen- 
tations labelled by J z (ii) CP with representations 
labelled by <?(+) and u(— ) and (iii) C3?. Here C 
interchanges Q and Q, P is parity and K is a ro- 
tation of 180° about the mid-point around the y 
axis. The C3? operation is only relevant to clas- 
sify states with J z — 0. The convention is to label 
states of J z = 0, 1, 2 by S, n, A respectively. The 
ground state (£+) will have J z = and CP = +. 

The exploration of the energy levels of other 
representations has a long history in lattice stud- 
ies [13114) . The first excited state is found to be 
the n u . This can be visualised as the symmetry 
of a string bowed out in the x direction minus the 
same deflection in the —x direction (plus another 



component of the two-dimensional representation 
with the transverse direction x replaced by y), 
corresponding to flux states from a lattice opera- 
tor which is the difference of U-shaped paths from 
quark to antiquark of the form n — U. 

Recent lattice studies ^S] have used an asym- 
metric space/time spacing which enables excited 
states to be determined comprehensively. These 
results confirm the finding that the H u excitation 
is the lowest lying and hence of most relevance to 
spectroscopy. 

From the potential corresponding to these ex- 
cited gluonic states, one can determine the spec- 
trum of hybrid quarkonia using the Schrodingcr 
equation in the Born-Oppenheimer approxima- 
tion. This approximation will be good if the 
heavy quarks move very little in the time it takes 
for the potential between them to become es- 
tablished. More quantitatively, we require that 
the potential energy of gluonic excitation is much 
larger than the typical energy of orbital or radial 
excitation. This is indeed the case ^31 > especially 
for b quarks. Another nice feature of this ap- 
proach is that the self energy of the static sources 
cancels in the energy difference between this hy- 
brid state and the QQ states. Thus the lattice 
approach gives directly the excitation energy of 
each gluonic excitation. 

The n„ symmetry state corresponds to exci- 
tations of the gluonic field in quarkonium called 
magnetic (with L PC = 1 H ) and pseudo-electric 
(with 1 h ) in contrast to the usual P-wave orbital 
excitation which has L PC = Thus we ex- 

pect different quantum number assignments from 
those of the gluonic ground state. Indeed com- 
bining with the heavy quark spins, we get a de- 
generate set of 8 states: 



L PC 


J PC 


J PC 


J PC 


J PC 


1-+ 


1— 


o-+ 


1-+ 


2-+ 


1— 


1++ 


o+- 


1+- 


2+~ 



Note that of these, J PC = 1"+, 0+ - and 2+~ 
are spin-exotic and hence will not mix with QQ 
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states. They thus form a very attractive goal for 
experimental searches for hybrid mesons. 

3.2. Hybrid decays 

Within this static quark framework, one can 
explore the decay mechanisms. One special fea- 
ture is that the symmetries of the quark and 
colour fields about the static quarks must be pre- 
served exactly in decay, hence the light quark- 
antiquark pair produced must respect these sym- 
metries. This has the consequence that the de- 
cay from a H u hybrid state to the open-6 mesons 
(BB, B*B, BB*, B*B*) will be forbidden PIT] 
if the light quarks in the B and B* mesons are in 
an S-wave relative to the heavy quark (since the 
final state will have the light quarks in either a 
triplet with the wrong CP or a singlet with the 
wrong J z ). The decay to B** -mesons with light 
quarks in a P-wave is allowed by symmetry but 
not energetically. 

In the heavy quark limit, the only allowed 
decays are when the hybrid state de-excites to 
a non-hybrid state with the emission of a light 
quark-antiquark pair. Since the H u hybrid state 
has the heavy quark-antiquark in a triplet P-wave 
state, the resulting non-hybrid state must also 
be in a triplet P-wave since the heavy quarks do 
not change their state in the limit of very heavy 
quarks. Thus the decay for b quarks will be to 
Xb + M where M is a light quark-antiquark me- 
son in a flavour singlet. This proceeds by a dis- 
connected light quark diagram and it would be 
expected ^7J that the scalar or pseudoscalar me- 
son channels are the most important (i.e. they 
have the largest relative OZI-rule violating con- 
tributions). This transition can be estimated on a 
lattice when the initial and final energies are sim- 
ilar. This is the case [7] for the II U de-excitation 
to ground state gluonic field plus /o meson when 
the interquark separation is around 0.2 fm which 
allows a lattice evaluation of the hadronic transi- 
tion strength - see fig. H3 

Note that in fig. © for R = 0.1 fm the II„ 
state is unstable to decay and, since we are us- 
ing Nf — 2 flavours of sea-quark, this decay is 
enabled and we might expect to see the lower en- 
ergy (i.e. £+ + /o) instead for this energy level. 
In practice, as we shall find later, the transition 
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Figure 6. The potential energy V(R) (in lattice 
units with a=0.0972 fm) versus quark separation 
R in fm for 2 flavours of sea quark. The energies 
are given [2] for the ground state and first excited 
gluonic state and for the two body state of ground 
state potential plus scalar meson (/o) in a P-wave 
with the minimum non-zero momentum. The on- 
shell transition can be evaluated when R ss 0.2 
fm. 



strength is quite small with x « 0.01 so the ad- 
mixture of the lighter state is very small indeed 
(it is of order x 2 ) and hence it will only be sig- 
nificant at very large i-values, much larger than 
those used here. 

Indeed the dominant transition is found to be 
with M as a scalar meson and the signal for the 
hadronic transition strength x is shown in fig. [7J 
Making use of the transition strength evaluated 
for a range of R- values around 0.2 fm and fold- 
ing in the wavefunction overlaps, we can estimate 
the decay width of the hybrid spin-exotic meson 
H (e.g. J PC = 1 h ) which has II K excited glue. 
This corresponds to a width of around 100 MeV 
for the transition H — > Xb + fo, whereas when 
M is an rj or rf meson the transition strength is 
found to be less than a few MeV. There will be 
modifications to this analysis coming from cor- 
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Figure 7. The transition matrix element xt for 
n„ — ► S+/o with momentum n = (1,0,0) versus 
t/a. Here R — 0.1 to 0.6 fm is represented by 
symbols: fancy square, diamond, +, octagon, x, 
square. The line represents a linear fit to the 
R = 0.2 fm case. 



rections to the heavy quark limit (of order 1 / ttiq 
where itiq is the heavy quark mass) which might 
allow hybrid meson transitions to BB, etc, but 
these have not been evaluated yet. 

In this heavy quark (or static) limit, the spin- 
exotic and non spin-exotic hybrid mesons are de- 
generate. For the latter, however, the interpre- 
tation of any observed states is less clear cut, 
since they could be conventional quark antiquark 
states. Moreover, the non spin-exotic hybrid 
mesons can mix directly (i.e. without emission 
of any meson M) with conventional quark anti- 
quark states once one takes into account correc- 
tions (of order 1/Mq) to the static approximation 
applicable for heavy quarks with physical masses. 

It is encouraging that the decay width comes 
out as relatively small, so that the spin-exotic hy- 
brid states should show up experimentally as suf- 
ficiently narrow resonances to be detectable. This 
decay analysis does not take into account heavy 
quark motion or spin-flip and these effects will 
be significantly more important for charm quarks 



than for 6-quarks. 

4. CONCLUSIONS 

It is possible to extract information about 
hadronic decays from euclidian lattice studies. 
The most rigorous method is to measure the small 
energy shifts in the two body states precisely as 
the lattice volume is varied. This then allows the 
phase shift to be extracted which gives informa- 
tion on resonance properties. This energy shift is 
only enabled in dynamical fermion simulations if 
quark- antiquark production is needed in the de- 
cay. This makes this a very difficult route to fol- 
low in practice. 

We have also sketched an alternative which 
is less rigorous but will enable estimates to be 
obtained more readily. This direct method was 
tested in the case of p — > tttt where it gave consis- 
tent results and, moreover, results in agreement 
with experiment. This is a good indication that 
the method is reliable in practice. The method 
works best when the lattice transition amplitude 
,x, is relatively small and it is encouraging that 
the method gives good indications when ax « 

0. 06. We then presented an application to hybrid 
meson decays: a case where the experimental re- 
sult is not known, but where the magnitude of 
x is indeed small. Our result is of experimental 
relevance: the knowledge of the magnitude of the 
decay width confirms that these hybrid states are 
likely to be accessible to straightforward study. 
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